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ABSTRACT 

Context. The isolated neutron star (INS) 2XMM J104608.7-594306 is one of the only two to be discovered through their thermal 
emission since the ROSAT era. Possibly a remnant of a former generation of massive stars in the Carina nehula, the exact nature of 
the source is unclear, and it might be unique amongst the several classes of galactic INSs. 

Aims. In a first dedicated XMM-Newton observation of the source, we found intriguing evidence of a very fast spin period of 
P ~ 18.6 ms at the 4cr confidence level. Moreover, spectral features in absorption have also been identified. We re-observed 
2XMM J104608.7-594306 with XMM-Newton to better characterise the spectral energy distribution of the source, confirm the candi¬ 
date spin period, and possibly constrain the pulsar spin-down. 

Methods. We used the two XMM-Newton observations of 2XMM J104608.7-594306 to perform detailed timing and spectral X-ray 
analysis. Both the spin-down rate and the energy of the spectral features provide estimates on the neutron star magnetic field, which 
are crucial for investigating the evolutionary state of the neutron star. 

Results. Statistically acceptable spectral fits and meaningful physical parameters for the source are only obtained when the residuals at 
energies 0.55 keV and 1.35 keV are taken into account by the spectral modelling. While the former can result from the inhomogeneous 
temperature distribution on the surface of the neutron star or can be related to a local overabundance of oxygen in the Carina nebula, 
the one at 1.35 keV is only satisfactorily accounted for by invoking a line in absorption. In this case, the best-fit neutron star atmosphere 
models constrain the hydrogen column density, the effective temperature, and the luminosity of the source within Ah = (2.5 - 3.3) x 
10^' cm“^, Teff = (6 - 10) X 10^ K, and Lx = (1.1 - 7.4) x 10^^ ergs“* *. The implied distance is consistent with a location in (or in 
front of) the Carina nebula, and radiation radii are compatible with emission originating on most of the surface. Non-thermal X-ray 
emission is ruled out at levels above 0.5% (3cr) of the source luminosity. Unfortunately, the second XMM-Newton observation proved 
inconclusive in terms of confirming (discarding) the fast candidate spin, providing an upper limit on the pulsed fraction of the source 
that is very close to the limiting sensitivity for detecting the modulation found previously. 

Conclusions. In the absence of an unambiguous period determination and an estimate of the magnetic field, the nature of the source 
remains open to interpretation. Its likely association with the Carina cluster and its overall spectral properties (only shared by a handful 
of other peculiar INSs) disfavour a standard evolutionary path or one in which the source was previously recycled by accretion in a 
binary system. The star 2XMM J104608.7-594306 may be similar to Calvera (IRXS J141256.0-I-792204), a neutron star for which 
the scenario of an evolved anti-magnetar has been discussed. A better age estimate and deeper radio and y-ray limits are required to 
further constrain the evolutionary state of the neutron star. 
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1. Introduction 

In the standard scenario of pulsar evolution, a neutron star is born 
rotating fast. As a consequence of rotational energy losses, it is 
expected to evolve towards longer periods until electron-positron 
pairs can no longer be produced; the radio emission ceases and 
the pulsar crosses the so-called death line. The radio emission is, 
however, energetically unimportant relative to the power emitted 
at high energies. While young, these objects are strong sources 

Send offprint requests to: A. M. Fires 

* Based on observations obtained with XMM-Newton, an ESA sci¬ 
ence mission with instruments and contributions directly funded by 
ESA Member States and NASA (Target 2XMM J104608.7-594306, ob- 
sids 0650840101, 0691970101). 


of X-ray and y-ray emission with radiative processes that include 
non-thermal emission from charged particles that were acceler¬ 
ated in the pulsar magnetosphere and thermal emission from the 
hot surface. However, after a maximum of ~ 10^ yr, pulsars have 
probably exhausted their internal source of particle creation and 
acceleration and cooled down to temperatures undetectable in 
X-rays. 

This relatively simple picture has been challenged by the dis¬ 
covery of peculiar groups of isolated neutron stars (INSs) that 
are radio quiet or transient. Their (mostly unpredicted) proper¬ 
ties signi ficantly differ from those d isplayed by standard pul¬ 
sars (see lKasMl2010t lHarding||2013L for an overview and ref¬ 
erences). Among them are magnetars (anomalous X-ray pulsars 
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(AXPs) and soft gamma repeaters (SGRs)), the central compact 
objects in supernova remnants (CCOs), the rotating radio tran¬ 
sients (RRATs), and the RO^Ar-discovered thermally-emitting 
INSs, also known as the Magnificent Seven (M7). 

In particular, the M7 constitute a nearby {d < 1 kpc) group 
of middle-aged (0.1-1 Myr) neutron stars, observable through 
their b right and purely thermal X-ray emission (see, for ex¬ 
ample, lHaberll[2007t lKaDlanl[2008t iTurollal[20091 for reviews). 
Amounting to abo ut half of all youn g INSs known in the so¬ 
lar neighbourhood (IPodov et al.ll2003l) . they could be as numer¬ 
ous as radio pulsars, with consequences for th e total number of 
neutr on stars populating the Milky W ay (see Keane & Kramei 


2008). Timing studies in X-rays (see iKanlan & van Kerkwiik 


201 11 and references therein) have shown that the M7 rotate 
more slowly (P ~ 3 - 11 s) and have higher magnetic field inten¬ 
sities, Bdip ~ (1 - 3) X 10'^ G, than the bulk of the (rotation- 
powered) radio pulsar population. Also at variance with pul¬ 
sars detected at high energies, the X-ray luminosity of the M7 
is in excess of the spin-down power. To a certain extent, such 
properties resembl e those of the you ng and energetic magnetars 
(see, for example, lMereghettill2008l for a review). On the other 
hand, the M7 experience less dramatic spin-down and do not 
show the complex phenomenology usually observed in magne¬ 
tars (e.g. emission of flares, bursts, quasi-periodic oscillations, 
spectral variability, timing noise, etc). 

Another group consisting of radio-quiet, t hermally emit¬ 
ting INSs, is that of CCO s (see, for example, 
iGotthelf & Haloed 120101 iHalnern & Gotthell 


de Lucal [20081 

20101 for re¬ 


views). CCOs were first observed as point-like, radio-quiet X- 
ray sources, located near the geometrical centre of supernova 
remnants. Of the dozen objects known at present, including other 
candidates, at least three fit in the anti-magnetar scenario, in 
which these sources are young (aged 10^ yr to 10^ yr) and weakly 
magnetised (Bdip = 10*° - 10** G) neutron stars that are experi¬ 
encing very low spin-down (IGotthelf et al.l20f3al and references 
therein). 

If the weak dipolar held of anti-magnetars remains constant, 
they are confined to a region of the P - P diagram of INSs that 
is devoid of other pulsars, between those occupied by the nor¬ 
mal and millisecond (recycled) part of the population. The fact 
that no old (orphaned, or without a supernova remnant) CCO is 
recognised in t his region of the P - P diagram in either radio or 
X -rav surveys dHalnern & Gotthell 2010l ^tthelf et ^l2013bl 


iLuo et alj2015t see also lBogdanov et al.l2014l who searched for 

CCO descendants among supposedly old radio pulsars in super¬ 
nova remnants) may favour the al ternative scenario of held burial 
by hypercritical accretio n (e.g. ICheyalieilll989t iGeppert et alJ 
Il999t IVigand et ai]|2012l) . In this scenario, the original neutron 
star magnetic held is submerged by fallback accretion of super- 
noya ejecta, remaining hidden for seyeral ~ 10^ yr to 10^ yr. As 
a result, the neutron star shows, in its early eyolution, the typ¬ 
ical anti-magnetar behayiour. As the held re-emerges, the pul¬ 
sar quickly joins the rest of the population, moying upwards in 
the P - P diagram; its subsequent eyolution would then follow 
according to the original field intensity. Residual thermal lumi¬ 
nosity in the first ~ 10“* yr, as well as enhanced cooling after 
~ 10^ yr, compared to its peers of similar characteristic ages, are 
expected to result from the accreted light-element enyelope (e.g. 
lYakoyley & Pethi^l2004l) . 

For seyeral years, considerable efforts ha ye been made to 
disco y er new thermally e m itting INSs (e . g. Chieregato et alJ 


[200l 

2XMM J104608.7-594306 (hereafter J1046), is one of only 
two to be discoyered through their thermal emission and 


Table 1. Instrumental configuration and duration of the EPIC 
scientific exposures of the XMM-Newton AO 11 observation of 


2XMM J104608.7-594306 


Instrument 

Start time 
(UTC) 

Mode 

Duration 

(s) 

pn 

2012-12-20T19:45:17 

SW 

87277 

MOSl 

2012-12-20T19:40:02 

FF 

87532 

MOS2 

2012-12-20T19:40:02 

FF 

87520 


Notes. The EPIC cameras were operated in imaging mode and the thin 
filter was used (obsid 0691970101). 


lack o f obyious counterparts since the ROSAT era (iPires et al.l 
l2009al) . Possibly a remnant of a former generation of mas- 
siye stars in the Carina nebula (iTownsley et al.[ [201 ll) . its ex¬ 
act nature is not clear, and it may be unique. The other 
source, IRXS J141256.0-1-792204, also known as Calyera 
([Rutledge et alJ200^ . is a relati yely bright R(9RA r source with a 
short spin period of P ~ 59 ms (iZane et alJl201 ih . Calyera is the 
only INS to date to be recognised as an eyolyed anti-magnetar, 
with a magnetic field of Bdip = 4.4 x 10** G tha t is, possibly, 
regaining its original strength (iHalpern et al.ll20l'^ . 

Thanks to its proximity (angular distance of 8.5') to the 
well-studied Eta Carinae system, J1046 was serendipitously de- 
tected on many o ccasions by the XMM-Newton Obseryatory 
(iJansen et al.l200l|) . The analysis of the serendipitous data shows 
no significant long-term yariability in either the X-ray flux or i n 
the oyerall spectral properties of the source ([Pires et al.l[2009bl) . 
Follow-up obseryations in the optical with the European South¬ 
ern Obseryatory Very Large Telescope (ESO-VLT) set a deep 
limit on the brightness of the optical counterpart, my > 27 (2cr), 
and a yery high X-ray-to-optical flux ratio of log(Px/Ev) > 3.8. 
At present, no radio or y-ray counterparts are known. 

As reyealed by a first dedicated XMM-Newton obseryation, 
performed on 2010 December 6 (A09, as reference for the text), 
the spectral energy distribution of J1046 is soft and purely ther¬ 
mal with sign ificant deyiations from a blackbody continuum of 
kT ~ 135 eV (iPires et al.l2012[) . While these properties are remi- 
niscent of those of the middle-aged M7, the XMM-Newton obser¬ 
yation also reyealed intriguing eyidence of a yery fast rotation, 
at P ~ 18.6ms. Such a rapid spin is difficult to reconcile with 
the spectral properties of the source. 

We re-obseryed J1046 with XMM-Newton two years later to 
further improye the characterisation of the source’s spectral en¬ 
ergy distribution, confirm the candidate spin period, and eyentu- 
ally constrain the pulsar spin-down. We report here the results of 
this obseryation in detail. We also reanalysed the A09 dataset to 
consistently compare our findings in a joined timing and spectral 
X-ray analysis. The paper’s outline is as follows: in Section[2we 
describe the new XMM-Newton obseryation and the data reduc¬ 
tion. Analysis and results are presented in detail in Section[3 We 
examine the eyolutionary state and possible nature of the neutron 
star in light of the obseryed properties of the known galactic INS 
population in Section |4| The main conclusions and results are 
summarised in Section[5] 

2. XMM-Newton observation and data reduction 

For details on the instrume ntation set-up and execution of the 
A09 obseryation, we refer to lPires et al.l(l201^ . The new XMM- 
Newton obseryation of J1046 (AOll) was carried out on 2012 
December 20 for a total exposure time of 87.714ks. Table [1] 
contains information on the scientific exposures and instrumen- 
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Table 2. Parameters of 2XMM J104608.7-594306, as extracted from the AOl 1 (and A09) XMM-Newton EPIC images 


Parameter 

pn 

MOSl 

MOS2 

EPIC 

EPIC (A09) 

Detection likelihood 

6596 

1613 

2280 

10567 

8952 

Counts 

5.54(9) X 10^ 

1.55(5) X 10^ 

1.71(5) X 10^ 

8.85(12) X 10^ 

7.83(11) X 10^ 

Counts (0.2-0.5keV) 

1.11(4) X 10^ 

2.70(22) X 102 

2.25(18) X 102 

1.61(5) X 10^ 

1.53(5) X 10^ 

Counts (0.5 - l.OkeV 

3.50(7) X 10^ 

8.3(4) X 102 

1.00(3) X 10^ 

5.36(9) X 10^ 

4.73(9) X 10^ 

Counts (1.0 - 2.0 keV) 

9.2(4) X 102 

4.44(29) X 102 

4.81(27) X 102 

1.86(5) X 10^ 

1.56(5) X 10^ 

Counts (2.0 - 4.5 keV) 

11 + 11 

0 + 2.0 

0.0 + 2.0 

13 + 11 

0 + 4 

Counts (4.5 - 12keV) 

2 + 6 

0 + 2.0 

4 + 4 

5 + 7 

12+10 

Rate (10-2 s-') 

9.18(15) 

2.23(8) 

2.18(6) 

13.67(18) 

13.27(19) 

HRi 

0.518 + 0.015 

0.51+0.03 

0.633 + 0.027 

0.539 + 0.012 

0.511+0.013 

HRz 

-0.585 + 0.015 

-0.31+0.04 

-0.35 + 0.03 

-0.485 + 0.012 

-0.505 + 0.014 

HR 3 

-0.976 + 0.023 

-1.000 + 0.007 

- 1.000 + 0.010 

-0.986 + 0.012 

- 1.000 + 0.006 

RA 

10 46 08.75(15) 

10 46 08.77(26) 

10 46 08.76(21) 

10 46 08.7(6)T 

10 46 08.8(4) 

DEC 

-59 43 05.32(15) 

-59 43 04.97(26) 

-59 43 05.07(21) 

-59 43 06.7(5)^ 

-59 43 06.8(5) 


Notes. Counts and rates are given in the total XMM-Newton energy band (0.2 - 12keV), unless otherwise specified. ^The marked EPIC source 
coordinates are astrometrically corrected, using as reference the 2MASS catalogue (see text). The corresponding Icr errors take the astrometric 
errors into account. 


tation set-up of the EPIC pn (IStriider et al.l 1200 111 and MOS 
dTurner et ai.ll200lll detectors. We adopted the thin filter for all 
EPIC instruments, given its better response at soft X-ray ener¬ 
gies and the very low (< 1 %) expected level of spectral pile- 
up at aimpoint. The pn camera was operated in small window 
(SW) mode, providing high time resolution (6 ms) at the ex¬ 
pense of a reduced exposure (deadtime of ~ 29%). Although 
the MOS cameras operating in timing uncompressed mode could 
provide sufficient resolution (1.75 ms) for investigating the pe¬ 
riodicity at P ~ 18.6 ms, calibration u ncertainties between 
the EPIC detectors (e.g. IZane et aT]l201 ih and the impossibil¬ 
ity of using the data for spectral analysis or for astrometrical 
corrections motivated us to choose the full-frame (EE) imag¬ 
ing mode instead. Standard data reduction was performed with 
SAS 13 (xmmsas_20131209_1901-13.5.0), using the latest cali¬ 
bration files. We processed the MOS and pn raw event files using 
the EPIC meta tasks emchain and epchain, respectively, applying 
default corrections. We ensured that the pn event file was clean 
of unrecognised time jumps (i.e. those uncorrected by standard 
SAS processing) by checking the output files of the epproc rou¬ 
tine. Additionally, we reprocessed and reduced the A09 obser¬ 
vation consistent with the new dataset. 

A few minor background flares were registered during the 
observation, effectively reducing the exposures in the MOS 1 and 
MOS2 cameras to 79 ks and 84 ks, respectively. In the pn camera, 
the background level for the duration of the observation was be¬ 
low what is recommended by standard filtering of ‘good-time in¬ 
tervals’. Therefore, the effective observing time (60 ks) was only 
reduced by the camera livetime in SW mode. 

Eor the analysis, we filtered the event lists to exclude bad 
CCD pixels and columns, as well as to retain the predefined pho¬ 
ton patterns with the highest quality energy calibration. Unless 
otherwise noted, single and double events were selected for pn 
(pattern < 4) and single, double, triple, and quadruple for MOS 
(pattern < 12). We defined the source centroid for each EPIC 
camera with the task eregionanalyse, which results at optimal 
radii of ~ 20" (0.3 - 2keV). Background circular regions of size 
60" to 100 " (depending on the camera) were defined off-source 
on the same CCD as the target whenever possible. 

The detected source count rates, hardness ratios, and other 
parameters based on a maximum likelihood PSP fitting are listed 
in Table|2] These are determined with the SAS task emidetect on 
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Fig. 1. Best-fit column density Ah as a function of blackbody temper¬ 
ature kT, for each individual (pn/MOSl/MOS2) spectrum used in the 
analysis (see legend). Also shown are the results by fitting all data si¬ 
multaneously, as well as results of the stacked (combined) spectrum. 
AOl 1 datasets are highlighted in red. Error bars are Icr confidence lev¬ 
els. 


images created for each camera and in each of the five prede¬ 
fined XMM-Newton energy bands. Eor comparison, we also list 
the source parameters as determined from the EPIC detection in 
the A09 observation. The errors are nominal Icr statistical un¬ 
certainties. There is no sign of variability in the source’s overall 
properties; typically, J1046 displays a very soft energy distribu¬ 
tion, with no counts above 2 keV, and count rates and hardness 
ratios are consistent within errors between the two XMM-Newton 
observations (see also Sect. EB- 

We used the task eposcorr to refine the astrometry, by cross- 
correlating the list of EPK T X-ray source positio ns with those 
of near-infrared (2MASS. ISkrutskie et al.ll2006l) objects lying 
within 15' from J1046. Based on a number of 60 X-ray/infrared 
matches, we found a slight offset in right ascension, not signif¬ 
icant within the astrometric errors, and an offset of 1.6" + 0.5" 
in declination. We double-checked the offsets with eposcorr us¬ 
ing other catalogues, in particular the Chan dra Carina Cluste r 
Project catalogue of X-ray sources (C CCP; BropseLaL 201 Ih . 
the Guide Star Catalog (GSC 2.3.2; iLasker et~^ 2008), and 
the Naval Observator y Merged Astrometric Dataset (NOMAD; 
IZacharias et al.ll2004l) . The EPIC source position (Table |B, up¬ 
dated accordingly, was found to be consistent within errors with 
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previous determinations. For consistency, we revised the source 
position in the A09 observation also by cross-correlating the 
EPIC source list with the 2MASS catalogue. We found offsets 
of -2.6" + 0.4" in right ascension and 1.6" + 0.5" in declination 
(46 matches). The updated source coordinates in the A09 epoch 
are also listed in Table |2] 

The statistics for the EPIC lightcurves (corrected for bad 
pixels, deadtime, exposure, as well as background counts, and 
binned into 875 s intervals) show the 3cr upper limits for the 
r.m.s. fractional variation of 0.12, 0.11 and 0.24 for pn, MOSl, 
and MOS2, respectively. On the basis of the same lightcurves, 
the reduced chi-square, assuming a constant flux, is 1.02 (pn, 99 
d.o.f.), 0.97, and 0.71 (MOSl and MOS2; 98 d.o.f.). 

3. Analysis and results 

3.1. Spectral analysis 

The analysis of the EPIC data is based on source and background 
spectra extracted from regions as described in Section|2] together 
with the respective response matrices and ancillary files created 
for each of the EPIC cameras. Owing to the frequent readout of 
pn in SW mode, detector noise dominates the low-energy count 
distribution. We therefore restricted the analysis of the pn camera 
to photons with energies above 0.35 keV. The energy band of the 
MOS cameras is restricted above 0.3 keV, which is in accordance 
with the guidelines and calibration status of the instrument. In¬ 
cluding the A09 observation, the available data amount to a total 
of six spectra and ~ 1.7 x 10"^ counts (0.3 - 12keV), of which 
10 % can be ascribed to background. Eor each spectrum, energy 
channels within 0.3/0.35 keV and 2keV were rebinned accord¬ 
ing to a minimum number of 25 counts per spectral bin. We also 
took care not to oversample the instrument energy resolution at 
a given bin by more than a factor of three, which is especially 
important at low energies; this is done by fixing the parameter 
oversample = i in the SAS task specgroup. __ 

To fit the spectra we used XSPEC 12.8.2 (lArnaudll996h . Un¬ 
less otherwise noted, the fit parameters were allowed to vary 
freely and within reasonable ranges. The p hotoelectric a bsorp- 
tion model and elemental abundances of IWilms et alJ (l2000i 
tbabs in XSPEC) were adopted to account for the interstellar ab¬ 
sorption. Owing to the intrinsically soft energy distribution of 
J1046 and its location in the galactic plane (b - -0.6°), as well 
as to uncertainties resulting from the choice of abundance table 
in XSPEC, we estimate that the amount of intervening material 
derived from the spectral fits is uncertain by 15%. 

Eirst, we fitted each of the six spectra individually, using a 
simple absorbed blackbody model. The resulting best-fit param¬ 
eters (Table[3]and Eig.[T]i argue against any significant variability 
between the two XMM-Newton observations. However, the Ah 
and kT derived from the MOS spectra, even though statistically 
consistent, are systematically higher and softer, by 30% and 5% 
respectively, than those derived from pn. We then proceeded with 
the analysis using two different approaches, both relying on the 
assumption that the spectral parameters of J1046 are constant be¬ 
tween the XMM-Newton observations: (i) we performed simulta¬ 
neous fits of all spectra, where we allowed for a renormalisation 
factor to account for cross-calibration uncertainties between the 
detectors, and {ii) using the task epicspeccombine, we merged 
all six spectra and corresponding background and response files 
into one single, stacked dataset. While both approaches give a 
similarly well-constrained parameter space, the exercise shows 
that the best-fit Ah and kT, derived from the second approach, 
are biased and inconsistent with those from individual fits. On 



0.5 1 2 

Energy (keV) 



Energy (keV) 

Fig. 2. Results of spectral fitting of source 2XMM J104608.7-594306. 
We show the two pn spectra and folded best-fit models, with residuals. 
Top: Simple absorbed blackbody. Bottom: Best-fit neutron star atmo¬ 
sphere aid - 2kpc (S = lO*^ G, - 10® K, Ah = 2.6 x lO^* cm-^). 
Two Gaussian lines in absorption are added to the model (with energies 
around 0.55 keV and 1.32 keV). 

the other hand, the results from the first approach, while natu¬ 
rally dominated by the better statistics of the two pn spectra, are 
consistent with most measurements. We henceforth adopted for 
the spectral analysis simultaneous fits, including all EPIC data 
from the two XMM-Newton observations accordingly. 

We summarise the results of the analysis in Tables 0] |5] 
and|6] Table 0] contains the results of spectral fits where we as¬ 
sumed a (single- or double-temperature) blackbody continuum, 
and tested for different element abundances and for the presence 
of absorption lines. The fitted models are labelled (A) to (E), as 
reference for the text (see caption). To test for overabundance of 
elements with transitions in the range of 0.3 keV to 2 keV (C, N, 
O, Ne, Mg, etc), we adopted the variable photoelectric absorp¬ 
tion model vphabs and restricted abundance values in the range 
1 < Z < 5 in solar units. 

The best-fit parameters of (more physically-motivated) neu¬ 
tron star atmosphere models are in Tables |5] and |6] where they 
are labelled (a) to (e). We tested magnetised and non-magnetised 
model atmospheres, consisting of fully- and partially-ionized 
hydrogen, as well as of carbon and other mid-Z e lements (the 
mode l s nsa, nsmaxg, and c arbatm in XSPEC: see IZavlin et al.l 
119961: iPavlov et al.l 1 19951: iHo et al.l 120081: iMori & HolliM 
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Table 3. Best-fit spectral parameters and count rates of 2XMM J104608.7-594306 per XMM-Newton observation and EPIC camera 


Ref. 

MJD 

(days) 

Ah (xlO^^ cm 

-") 


kT (eV) 


Net count rate (xlO "^s ^) 

pn 

MOSl 

MOS2 

pn 


MOSl 

MOS2 

pn 

MOSl 

MOS2 

A09 

55536.531753 

9 7Q + U . Z4 

^•^^- 0.23 

9 ' 7 + 0.8 

^‘- 0.7 

9 O + O . / 

135(3) 


128(8) 

131(7) 

7.39(12) 

1.54(6) 

1 .66(6) 

AOll 

56282.326597 

9 91 + 0.27 
■^•^^- 0.26 

9 7 + 0.6 

4.3(5) 

133:^4 


131(6) 

125(5) 

7.39(12) 

1.63(5) 

1.73(5) 

Simultaneous (per camera) 

2.98(18) 

9 7 + 0.5 

4.1(4) 

134.2(2.4) 

130(5) 

127(4) 

7.39(17) 

1.58(8) 

1.70(8) 

Simultaneous (EPIC) 

3.20(15) 



133.1(1.9) 



21.34(20) 


Combined 

2.69(13) 



136.7(1.8) 



21.94(19) 



Notes. Errors are Icr confidence levels. The model fitted to the data is a simple absorbed blackbody (tbabs*bbody). Reducedvalues are within 
0.9 and 3.4, for 23 to 30 degrees of freedom (individual fits). Count rates in energy band 0.3 - 12keV are subtracted for the background. 

Table 4. Results of spectral fitting of blackbody models 


Model 

Tv (d.o.f.) 

Ah 

(xlO^' cm“^) 

(eV) 

kT^ 

(eV) 

G 

(keV) 

G 

(keV) 

(ergs^' cm“^) 

R“T 

(km) 

R“® 

(km) 

(A) 

1.98(169) 

3.20(15) 

- 

133.1(1.9) 

- 

- 

9.3!; “ X 

- 

3-8!;;:^ 

(B) 

1.34(167) 

5.3(3) 

28.6!“ 

117.8(2.2) 

- 

- 

6!“ X 10-'° 

4!4x 10®“ 

6.8!°;® 

(C) 

1.48(168) 

3.75(17) 

- 

121.3(2.1) 

- 

- 

1.5!°°x 10-'2 

- 

5-8!!, 

(D) 

1.08(165) 

9 17+0.22 
-0.20 

- 

157(5) 

- 

1 9QQ+0.018 

4.6!®4X 10-'® 

- 

^■^-0.4 

(E) 

1.03(165) 

2.16(20) 

- 

144(3) 

0 577+0-013 
' -0.012 


5.4'f X 10-'® 

- 

2.5(4) 

(F) 

1.02(165) 

T . 4+0.4 
-0.3 

36(4) 

150(5) 


1 97A+O.OI8 
-0.019 


3JQ+320A 

-1.30 

7 4+0.4 
^■^-0.3 


Notes. Errors are Icr confidence levels. All tbabs models assume I Wilms et^ ( |2000|) abundances. Models: (A) tbabs(bbody); 
(B) tbabs(bbody+bbody); (C) vphabs(bbody), with best-fit oxygen overabundance of Zq = 1.60(7) in solar units; (D) vphabs(bbody+gauss), 
with Zq = 1.64(12), best-fit Gaussian cr ~ 0.2keV, and line equivalent width of EW ~ 220 eV; (E) tbabs(bbody+gauss+gauss), with fixed Gaus¬ 
sian CTi =0.1 keV and cr 2 = 0.2 keV, and EWi ~ 90 eV and E'W 2 ~ 120 eV; (F) tbabs(bbody+bbody-(-gauss), with fixed Gaussian a 2 = 0.2 keV 
and E'W 2 ~ 210eV. The unabsorbed flux is in energy band 0.2 - 12keV. The radiation radii are estimated normalising the distance to the 
source to that of Eta Carinae, dear = 2.3 kpc. In order to constrain R" = 20 km, the source must be at 10 pc (R” 2km). In order to 

constrain R” = 20 km, the source must be at 150 pc (R” ~ 2 km). 


ISuleimanov et alJ l2014ll . Wherever appropriate, the magnetic 
field was held fixed at values B - 0, 10*^, 10*^ G (nsa, Table|5l); 
for nsmaxg, we tested 19 models with B - (0.01 - 30) x 10'^ G 
(only the results with < 2 are listed in Table |6]l. We first 
assumed canonical values for the neutron star mass and radius, 
Mns = 1.4Mo and Rns = 10 km; these parameters were also 
allowed to vary to check for an improved fit. The size of the 
emission region, with respect to the neutron star physical radius, 
can also be parametrised for nsmaxg models; we restricted val¬ 
ues to have Rem ^ Rns- Finally, we also added a second thermal 
component, varied the elemental abundances, and added lines 
to the continuum to obtain statistically acceptable fits and de¬ 
rive plausible physical solutions for the neutron star. For double¬ 
temperature nsa and nsmaxg models, we fixed the mass of the 
neutron star at 1.4Mo, let the radii (Rns or Rem) of the two com¬ 
ponents vary independently of each other, and limited the dis¬ 
tance to the source to within 10 pc and 50 kpc. The mass, mag¬ 
netic field, and distance normalisation of the second thermal 
component were set equal to those of the first. 

The analysis confirms the results of iPires et al.l (l2012h . with 
much improved count statistics. Overall, the spectrum of the 
source is soft and purely thermal. However, significant devi¬ 
ations from a Planckian shape are identified around energies 
0.55 keV and 1.35 keV (see top plot of Fig. |2]i. The residuals 
are present independently of the choice of the thermal contin¬ 
uum, abundance table, source or background extraction regions, 
as well as regardless of the EPIC instrument and XMM-Newton 
observation. 

The best-fit absorbed blackbody (A) has kT - (133.1 + 
1.9)eV and Ah = (3.20 + 0.15) x lO^^cm-^, with xl =1.981 for 
169 degrees of freedom (d.o.f.) (Table|4]l. The column density is 


consistent with the one to wards Eta Carina, Ah ~ 3 x 10^' cm“^ 
(iLeutenegger et al] l2003h . Assuming t he neutron st ar is at a 
comparable distance, dcnr - 2.3 kpc jSmi tiil l2006h . the red- 
shifted radiation radius derived from the blackbody fit is Roo = 
(3.8((Q 4 )(t//£/car) ' km. We note that the redshifted radiation radii 
of the M7, as derived from X-ray blackbody fits and distance es- 
tim ates, are in the range of roughly 2 km to 7 km (see, for exam¬ 
ple, [Kapi^|2008l). The residuals mentioned above mostly con¬ 
tribute to the high value of chi-square. 

Similarly, if the residuals are not taken into account by 
adding more spectral complexity, atmosphere models do not pro¬ 
vide good fits or plausible physical parameters for J1046. This is 
illustrated by fits (a) in Table|5]and|6] which consist of simply ab¬ 
sorbed nsa and nsmaxg models respectively, with fixed element 
abundances (no carbatm model provided a fit with < 10). 
The fits lead to inconsistent column densities, sizes of the emis¬ 
sion region, and distances; the quality is poor (xl - 1.6 - 1.9, 
for 166 to 169 d.o.f.) and several parameters are unconstrained. 
The effective temperature of magnetised nsa models is usually 
above 10® K as they consist of fully-ionized hydrogen atmo¬ 
spheres. The only nsmaxg models with < '2- consist of cold, 
Tes ~ (4 - 5) X 10® K, partially-ionized hydrogen atmospheres, 
under moderately strong magnetic fields of R ~ (2 - 3) x lO'® G. 
The large column density. Ah ~ (5 - 7) x 10^' cm“^, is incon¬ 
sistent with a best-fit distance of only ~ 60pc. We verified that 
other element atmospheres (C, O, and Ne) in nsmaxg do not im¬ 
prove the spectral fitting. 

To test a multi-temperature spectral distribution, we next 
tried adding one extra thermal component to the absorbed black¬ 
body and atmosphere models. With respect to (A), we obtain an 
improved fit quality (by A;^fy ~ 0.6), and blackbody temperatures 
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Table 5. Results of spectral fitting of neutron star atmosphere models (nsa) 


B xi (d.o.f.) 

A^h 

Teff,! 

7’eff,2 

ei EWi 

£2 

£ 1^2 M„, 

DOO 

^em,l 

poo 

^em,2 

log(Lx) 

d 

(G) 

(102'cm-2) 

(10® K) 

(10® K) 

(keV) 

(keV) 

(Mo) 

(km) 

(km) 

(ergs-') 

(kpc) 


(a) 

0 

1.92(169) 

4.85(11) 

4 38+^-^^ 

- 

- 

(b) 

lel2 

lel3 

0 

1.63(169) 

1.81(169) 

2.19(166) 

7 cq+0.13 

37o;8:!? 

^-0.23 

4 7 Q+O.IO 

2 84;8:i§ 

A 71 + 0.17 
^■'*-0.26 

30.0*4 

21.2!^; 

3.00!°^ 

10 

- 


lel2 

lel3 

1.68(166) 

1.96(166) 

^■^-1.7 

9 1 + 0.8 
■^■^-0.5 

24.58*°f 

13-52!™ 

: 

(c) 

0 

1.57(165) 

4.98(11) 

4 41 +0.26 
^*^^- 0.22 
33!| 

4 8 *' ® 

^■°-0.4 

- 

- 


lel 2 

1.32(165) 

9 70 + 0.13 
^■'^-0.20 

6 34-^0.19 

- 

- 


lel3 

1.89(165) 

- 

- 

(d) 

0 

1.39(166) 

2.80(3) 


15.1*f 

- 


lel 2 

1.15(166) 

1 40 + 0.22 

6(4) 

38!° 

22 

- 


lel3 

1.28(166) 

4 1*4 

^■■^-0.3 

10 

94 ? 8 :t 

- 

(e) 

0 

1.08(165) 

2.9(3) 

- 

0.548!g;g[’ 

0.549!;; 

0.545i°l| 


lel 2 

1.05(165) 

9 cq+ 0.21 
z..^y 03 

- 


lel3 

1.12(165) 

9 Qi +5.2o 

- 


1.4 
2 

2.5 
1.4* 
1.4* 
1.4* 

1.5 
1.9 
0.8 


13 

18 

25 

99 c+ 0.6 
^^•'^-0.9 

23 

23 

13 

18 

11 



1.346(15) 

0.3 

1.4* 

23 

^*^^^- 0.016 


1.209*° 

1.309*° 

012 

0.25 

0.2 

1.4* 

1.4* 

23 

23 

6 

6 

0.07 

1.327*° 


0.26 

1.4* 

13 


0.07 

1.320*° 


0.22 

1.4* 

13 


0.08 

1.352*° 

iJit 

.016 

0.28 

1.4* 

13 

- 


31.42 

34.44*1 

34.00 

32.34; 

33.82' 


8:lf 

Ue 


- 0.04 
704 + 0.04 
JZ. /77_q q24 
7 1 7 0+0.10 

-Q.Q9 

34 53 *° '^ 

71 ^ 

33 1 *°-^ 
35 ( 5 ) 
33 . 68 * 00 ^ 
32 05 *° '"^ 

32 . 87*^3 


32.75^ 


0.00 


0.17 

8 - 0 !^*^ 

2.58(4) 

n 99 C+O.OO 4 

70^8 

^•'^^^-0.005 

0.12 

7!^ 
0.10 
1 Q 9 + 0.25 
^■^*^-0.23 

13.8!°® 

7 70 + 0.09 
-’■■’'^-0.04 

0.654(6) 

1.997(19) 

1 'S74+0015 


Notes. Errors are Icr confidence levels (unconstrained parameters are shown with no errors for illustrative purposes). All tbabs models assume Iwilms et ^ (l2000h abundances. All listed models 
consist of hydrogen atmospheres. Models: (a) tbabs{nsa); (b) tbabs(nsa+nsa); (c) vphabs{nsa), with best-fit oxygen overabundance of Zq ~ 1.3 - 1.7 in solar units; (d) tbabs(nsa+nsa+gauss), 
with fixed Gaussian <T 2 = 0.2 keV; (e) tbabs(nsa+gauss+gauss), with fixed Gaussian <ti =0.1 keV and <T 2 = 0.2 keV. **' Parameter held fixed during fitting. 
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Table 6. Results of spectral fitting of neutron star atmosphere models (nsmaxg) 


B 

^!(d.o.f.) A^h 

Teff,! 

7’eff,2 

ei EWi 

ei 

EW 2 

7?” , 

em,l 

em,2 

log(Ax) 

d 

(G) 

( 10 ^’ cm“^) 

(10^ K) 

(lO^K) 

(keV) 

(keV) 

(Mo) 

(km) 

(km) 

(ergs^') 

(kpc) 


(a) 

2el3 

1.64(166) 

n 14 +U.U 3 

5 92^8:88 
-Q-Q5 

3.98 


3el3 

1.83(168) 

5.00 

(b) 

3el3 

1.86(166) 

-^•*^^-0.09 


(c) 

lelO 

1.69(167) 

5.46(21) 

3.88 


7el2 

1.92(167) 

7.03!Sf 

5.13 


lel3 

1.89(167) 

6 77 + 0-20 
-0.29 

5.40 


3el3 

1.56(167) 

5.15(5) 


(d) 

lelO 

lel 2 

1.34(164) 

1.70(164) 

2.80!; 

8.59 


1.3el2 

1.58(164) 

0 *79+0.13 
-0.27 

3.78!°-" 

3 

^- 0.15 

8.36 


2 el 2 

1.58(164) 



4el2 

1.47(164) 



7el2 

1.42(164) 

3.95(15) 



(e) 


lel3 

2el3 

3el3 

lelO 

lel 2 

1.3el2 

2 el 2 

4el2 

7el2 

lel3 

2el3 

3el3 


1.41(164) 

1.37(164) 

1.40(164) 

1.05(164) 

1.21(164) 

1.20(164) 

1.19(164) 

1.16(164) 

1.15(164) 

1.14(164) 

1.13(164) 

1.08(164) 


4 02^° '’' 

0.15 


4.15!; 

3.90!| 

2.49!| 

3.29!| 

3.21!| 

3.17!! 

3.05!; 

3.04!; 

3.03!; 

2.99!; 

3.21 


28 


io 

0.3 


-4 

+ 0.05 
-5 

3.2 
+ 0.16 

- 0.14 
+ 0.20 

■2_3 
3 + 0.16 
’-OJ 

3 

3 + 0.13 
-"-O.!! 

9 
3 
3 
3 
3 
3 

n 7 + 0.6 
- 0.5 

7.4(5) 
7.9(5) 
8 4^° "^ 
9278:^ 

9.4(5) 

9.5(5) 

9.6(5) 

8 A+O.V 


79.4 


1.1 

1.4 

1.4* 

1.4* 

1.4* 

1.4* 

1.4* 


- 0.537* 


113 

0 539*““'^ 

1.010 
1.015 
1.012 
0.550(11) 


0.541!; 

0.545!8 


0.08 

0.08 

0.08 

0.07 

0.08 


1.357! 

1.353*' 

1.352! 

1.345* 

1.313 


0.017 


0.29 

0.28 

0.28 

0.27 

0.20 


1.4* 

1.4* 

1.4* 

1.4* 

1.4* 


23 

13 

13 

13 

13 

13 

13 


10 

10 

- 


1 360^®-®^^ 

1 *?Q 1 + 0.008 
^■'^^^_n.no9 

1 369+0™ 

i.JSZ-ooii 

0.3 

0.4 

1.4* 

1.4* 

7.7(2.8) 

«c 2 

10 

- 


0.3 

1.4* 

13(4) 

10 

- 


0.3 

1.4* 

1 j 4 + 0.6 

8 . 6 !?-f 

10 

- 


1.382(10) 

0.3 

1.4* 

10 

- 


1 aco+O.Oll 
i.jjj_Ooio 

0.3 

1.4* 

r-- 

i? 

00 

00 

(N 

10 

- 


1.350(11) 

0.3 

1.4* 

10 . 0 *“-° 

10 

- 


1 346*°-°'^ 

1 999+0.015 
i.2)J2)_0.013 

0.3 

1.4* 

11.7(6) 

10 

- 


0.28 

1.4* 

12 4*°-^ 

- 

0.554(14) 

0.07 

1.349(16) 

0.27 

1.4* 

9 

_ 

0.525*°-°" 

0.08 

1.368*°-°" 

0.3 

1.4* 

11 

- 

n coo+ 8 :ffi 
^-^^o_q.009 
0.532+^!!;3 

0.08 

0.08 

1 366^8:8!t 

i.JUU_QQ 12 

1.368*°-°!! 

0.3 

0.3 

1.4* 

1.4* 

11 

10 


9 

10 
9 

8.5 

12 


31.92 

31.65 


-22 

+ 0.02 

-23 


« 1 

30.85 

_ 

31.21(7) 

_ 

31.69*°-® 

- 

31-78!°;° 

30.9!|° 

< 1.8 

99 1 + 0.8 

■2^-^-0.5 

2.6 

31(24) 

< 2.0 

32.54!°-® 

<2.3 

32.57*1® 


0 94^+0 010 
u.JU‘+_QQ 22 

0.254(16) 
0.446!°°f 
0 99 + 0.10 
'*■^^-0.16 


30 . 52!|'3 

30.87!“4‘'^ 

30.65!“2" 

30.82!°;® 

31-02!:S 

32 06*° '^ 
32.20!°" 
32.30*°-'° 
32.37_ 

9Z,.J)J_p .. 

99 4q+0 0» 
-QOS 
32 47+° ® 

-’^■^'-0.09 
99 41 +° *>^ 

32.51- 


m 


0.08 


0.06 

0.06 

0 051*°-°''' 

u.uji Q023 

0.09 

0.06 

0.08 

0.06 

1 . 001 ( 8 ) 

0.310(3) 

0.975(8) 

0.998(3) 

0.03 

0.05 

0.03 

0.06 

0.15 

0.844(7) 

0.579(5) 

0.702(6) 

0 7 Q " i + 0.007 
^*^-0.006 
0.855(7) 

1 rvoo+O.OOQ 

- 0.008 

0.996(8) 

0.956(8) 


Notes. Errors are Icr confidence levels (unconstrained parameters are shown with no errors for illustrative purposes). All tbabs models assume Iwilms et ^ ll2000h abundances. Only results of 
models with;^'J < 2 are shown. All listed models consist of hydrogen atmospheres. Models: (a) tbabs(nsmaxg); (b) tbabs{nsmaxg+nsmaxg); (c) vphabs{nsmaxg), with Zq ~ (1.2 - 1.3); 
(d) tbabs{nsmaxg+nsmaxg+gauss), with fixed Gaussian 0-2 = 0.2 keV; (e) tbabs{nsmaxg+gauss+gauss), with fixed Gaussian a-, =0.1 keV and 0-2 = 0.2 keV. **' Parameter held fixed during 
fitting. 
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of kT^ = 28.6:"2 ]eV and kT^ = (117.8+2.2)eV (B). However, 
the large column density of A^h = (5.3+0.3)xl0^'cm“^requires 
the model normalisation of the soft component to be large (the 
unabsorbed flux is nearly three orders of magnitude higher than 
that of a single-temperature blackbody). At 2.3 kpc, the radiation 
radius of the soft component, ~ 4,000 km, is inconsistent 
with that of a compact object. Even if the neutron star is in the 
foreground of the Carina complex, it should be at a maximum 
distance of 10 pc from the Sun to constrain < 20 km. Con¬ 
sidering the large Ah of the model, the low galactic latitude of 
the source, and a chance alignment with the Carina nebula, the 
solution is unlikely. 

A second thermal component in the atmosphere models 
worsens the chi-square (b). Nonetheless, in spite of a rather poor 
fit ~ 1.7 for 166 d.o.f.), one solution is more physical, the 
nsa model for B = lO'^G, with Ah ~ 2.84^°'2“ x 10^' crn^^, 
kT^ = 34+21 gY ^ 7^00 ^ i88.6!|®gv, at a best-fit distance 
of ~ 3.9 kpc. The best-fit radiation radii are, however, pegged 
at their boundary values of = 20 km and = 5 km. All 

other (b) models in Tables |5] and |6] show large high Ah, and 
distances below 500pc. 

In (C)/(c), we found that the agreement between model and 
data is improved below 0.7 keV (typically by ^ o.3 - 0.5), 
by allowing oxygen to be extrasolar, with Zq ~ 1.2 - 1.6 in 
solar units, depending on the choice of continuum. In general, 
there is no significant change in the best-fit parameters (the ex¬ 
ception is for the nsa model with for B - 10 ^^ g) Despite the 
improvement, all nsmaxg models (c) still show large Ah and dis¬ 
tances contradictorily below 100 pc, regardless of the magnetic 
field strength. The abundance of oxygen was the only one tightly 
constrained by the vphabs fits, which were found to be either in¬ 
sensitive to other elements (e.g. C, Ne, Si, or Fe), or to produce 
large overabundances (e.g. N, Mg, and Ca). In particular, ele¬ 
ments with lines around 1.35 keV were verified to be arbitrarily 
overabundant in the model. A variable absorption model, com¬ 
bined with a double-temperature blackbody, has no effect on the 
elemental abundances, which are insensitive to the fit. This is 
likely due to the fact that the two blackbody components cross 
around an energy of 0.5 keV, where the vphabs model tries to 
accomodate for the extrasolar amount of oxygen. 

Independently of the continuum - single or double¬ 
temperature, blackbody or atmosphere, with fixed or overabun¬ 
dant elements - we verified that the introduction of a Gaussian 
line in absorption around an energy of 1.35 keV significantly im¬ 
proves the quality of the spectral fitting. With respect to (C), 
the presence of the line improves the fit by A;^f2 o.4, result¬ 

ing in a less absorbed model, with Ah = 2.17+g2Q x 10^^ cm“2 
and kT — (157 + 5)eV (D). The oxygen abundance is simi¬ 
lar, Zq — (1.64 + 0.12), and a best-fit line energy is centred at 
e - 1 . 289 +gg[ 2 keV, with equivalent width EW ~ 220eV. Simi¬ 
lar parameters and chi-square values are obtained when, alterna¬ 
tively, the abundances are fixed, and two lines in absorption are 
added to the blackbody continuum (E). In this case, the best-fit 
line energies are centred at 0.577+QQ}2keV, and 1.361+gg[gkeV 
and have equivalent widths of 90 eV and 120eV. (To better con¬ 
strain the model parameters, the Gaussian cr of the two compo¬ 
nents were held fixed to 0.1 keV and 0.2 keV, respectively.) 

As for (B), the presence of a line at e = 1.270+QQ|gkeV re¬ 
sults in a lower flux normalisation and a column density con¬ 
sistent with that towards Carina, Ah = 3.4+g2 x lO^' cm“2 (p)- 
the fit quality is improved by A;^f2 .^0.3 and the temperature of 
the two components are higher (kr“ = 36 + 4 eV and kT^ - 
150 + 5 eV). Still, the neutron star should be at a maximum dis- 


Table 7. Best-fit parameters per observation and EPIC camera 



Nh''' 

TeS 

e 

d 


(Xl02') 

(xio^ K) 

(keV) 

(kpc) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

7 q+u.s 

2 3+1.2 

4 9+1.5 

-Q-5 

3 1+0-9 

-Q-S 

0 0+1.3 

Q ^ + 1.0 

^■*^-0.9 

9-i:li 

8.6!j;2 

9-4!!g 

9 1+1:8 

()36W™2Sj32(3) 

0.516+°°25 1 40(13) 

0.80(5), 1.24+°;"^ 

0 S44+0022 1 oqc+0.024 
^■‘^"+"+_0.019’ ^•^^‘^-0.025 

0.59+“°^.33(5) 
0.52!™, 1.32(5) 

1 5+u.f> 

^•'*’-0.5 

1 5+'-^ 

09;8:t 

9 7+0.6 

-Q -7 

1 4+0.9 

-^-0.5 


Notes. Errors are Icr confidence levels. The fitted model is (e) for a 
nsa with B = 10*2 G (see Tabled- The fitted spectra are (1) pn A09, 
(2) MOSl A09, (3) MOS2 A09, (4) pn AOll, (5) MOSl AOll, (6) 
MOS2 AOl 1. The reduced chi-square values range between 0.6 to 1.5. 
The column density is in cm“2 


tance of 150 pc to constrain < 20 km. Regarding the double¬ 
temperature atmosphere models, the presence of the absorption 
line at ~ 1.35 keV also results in solutions with more consistent 
physical parameters (d). In particular for nsa with B - OG and 
B - 10*2 G, the column density is Ah = (2.6 - 2.8) x 10^* cm “2 
and the distance is within d ~ 1.9kpc to 3.3 kpc, i.e. consis¬ 
tent with a neutron star in Carina. The best-fit solutions have 
kT” ~ 30 eV and kT^ = 120 - 170eV, with ~ 20 km and 

~ 5 km. 

Finally, we tested a combination similar to that in (E), i.e. a 
single-temperature model with two lines, for the nsa and nsmaxg 
models. These are models (e) in Table |5] and respectively, 
which are also the ones with the lowest-;^^^ values in the anal¬ 
ysis. Interestingly, all relevant parameters show best-fit values 
in a much narrower range. The energy of the lines are centred 
at e\ ~ 0.55 keV and 62 ~ 1.31 - 1.37keV, and the equiva¬ 
lent widths are EW\ ~ 72-81 eV and EW2 ~ 200 - 300 eV. 
As for (E), the Gaussian cr of the two components was fixed. 
Canonical Mns and values typically provided the least chi- 
square fits. The size of the redshifted emitting region is within 
8.5km and 13 km - overall consistent with the emission origi¬ 
nating from most of the surface. The magnetic field intensity is 
not constrained by the atmosphere models, as the least-j^f^ fits 
have B - (0 - 3) x lO'^G. In the lower plot of Figure |2] we 
show the two pn spectra folded with the best-fit nsa model (e), 
forB = 10*2 G. 

In Table|2]we summarise the best-fit parameters of this same 
model for the six individual EPIC spectra. The results are consis¬ 
tent overall. The most deviant parameters result from fitting the 
MOS2 spectrum of the A09 observa tion. Indeed, this spectrum 
was not considered in the analysis of iPires et al.l (l2012l) . 

The inclusion of power-law tails extending towards higher 
energies has no effect on the parameters of the least-;if 2 models. 
Considering again the best-fit nsa model (e), we found that non- 
thermal power-law componen ts with typical phot on indices of 
E = 1.7-2.1 (see, for example. lDe Luca et al.l2005l) contribute at 
most 0.5 % (3cr confidence level) to the luminosity of the source. 

3.2. Timing analysis 

To reassess the significance of the periodic signal detected at 
P ~ 18.6 ms, we restricted the timing analysis to the pn data and 
completely reanalysed the A09 dataset consistent with the new 
XMM-Newton observation. We considered events with photon 
patterns 12 or lower and tested different energy bands and radii 
of the source extraction region; the purpose was to unveil po- 
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Frequency (FIz) 


Fig. 3. Results of Z\ tests in the A09 (top) and AOl 1 (bottom) pn time 
series. The frequency range is Av ~ 87.7 Hz. Dashed horizontal lines 
show confidence levels of Icr to 4(T (68% to 99.994%) for the detection 
of a periodic signal, taking the whole frequency range of the search into 
consideration (see text). The frequency of the highest peak resulting 
from each search is labelled. Top: the energy band is 0.36 - 2.25 keV, 
the size of the extraction region is 18.6", and the number of photons is 
5188. Bottom', the energy band is 0.37-1.9 keV, the size of the extraction 
region is 19.7", and the number of photons is 5244. 


tential pulsations from the neutron star by minimising the back¬ 
ground contribution. The spatial and energy selection cuts pro¬ 
vide roughly Aiph ~ (4 - 6) x 10^ pn events in each observation. 

Each pn time series was Fourier-analysed in frequency do¬ 
main to search for the presence of pe riodic signals. The 
(Rayleigh) statistic dBuccheri et al.ll9^ was computed directly 
on the photon times-of-arrival; these were converted from the 
local satellite to the solar system barycentric frame using the 
SAS task barycen and the source coordinates in each pn obser¬ 
vation (Table |2|i. We adopted a step in frequency of 2y(Hz (or 
an oversampling factor of the expected width of the peaks of 
nearly 10), warranting that a periodic signal is not missed in our 
searches. The number of statistically independent trials in each 
Z^ test is estimated N - Avfexp ~ 8 x 10®, with Av ~ 87.7 Hz 
and fexp ~ 87 ks to 90 ks. In the absence of a periodic component, 
the mean and standard deviation of the Z^ values, with n - are 
equal to 2, as the Zj statistic is distributed as a with two de¬ 
grees of freedom. Therefore, the probability of obtaining a noise 
peak of a given power Z^ in N independent trials is given by 
a - A/'exp(-Z^/2); the confidence level of a detection can then 
be established as C = (1 - a) x 100%. 

The analysis of the A09 dataset yields only one statisti¬ 
cally significant periodic signal, at v, = 53.645804(3)Hz (P» = 
18.6407869(12)ms; see top plot of Fig. O. The power of the Zj 
statistic at v, is energy and signal-to-noise {S IN) dependent, be¬ 
ing measured at best (or most significantly) when photons with 
energy below ~ 0.35 keV (where the noise of the pn camera in 
SW mode dominates) are excluded from the analysis, or when 
the source extraction region is around an optimal value between 



Phase 


Fig. 4. A09 pn light curve, folded at R ~ 18.6 ms (the epoch of phase 
zero is defined at Tq = MJD 5536.00884 days). The energy band is 
0.36 - 2.25 keV and the size of the extraction region is 18.6". Two cy¬ 
cles are shown for clarity. The best-fit sinusoidal function is also shown 
(dashed line). 


17" and 22". These reflect the intervals where the S /N ratio is 
optimal, indicating an origin on the neutron star itself. 

The measured fluctuation of the power of the Z^ statistic at 
v,, Z^(v,) ~ 36 - 50, is consistent with that expected from a si¬ 
nusoidal modulation of amp litude pf - (12.5 + 1.9)% (see, for 
example. iPavlov et al.lll999l) . Indeed, this is the amplitude de¬ 
rived from a sinusoidal fit of the light curve of the pn camera, 
folded at P = P, (Fig. IH; the pulsed fraction is also consis¬ 
tent with the non-parametric determination of the signal strength, 
Pf - (13.1 + 1.6)%, computed directly from t he sample phases 
accord ing to the bootstrap method proposed bv ISwanepoel et aP 

dl^ . 

The inclusion of higher harmonics n > 2 in the Z^ test was 
found to be statistically insignificant. Extensive searches showed 
no other significant signal (above the Icr confidence level, or 
with Z^ > 34) across the frequency range of the analysis, and 
safety checks did not suggest an unknown inst rumental or back- 
ground origin for the periodic signal at v» (see lPires et al.ll201^ 
for details). We also verified that the periodic signal is always 
present at the same frequency, independently of the details on 
the processing of the raw event file (e.g. included calibration 
files, SAS version, randomisation in energy within a PI chan¬ 
nel, event filtering, or randomisation in time within the sampling 
detector time). 

In the AOl 1 dataset, however, no obviously significant can¬ 
didate frequency results from our Z^ tests. Extensive tests posi¬ 
tively exclude modulations with pf > 14% (Zj >51, 4cr) in the 
v = (0.01 - 87.73) Hz frequency range of the search. However, 
at the brightness level of J1046, a modulation with pf = 11% to 
13% may be missed in searches with too large a number of tri¬ 
als. As the background level in the second observation is higher 
and the exposure time slightly shorter, the pulsed signal at P, ms 
may have been more severely affected (or smeared out) by noise. 
Alternatively, the A09 signal is spurious and/or unrelated with 
the spin period of the neutron star (see discussion in Sect. 14. 11 1. 

By considering Av = 87.7Hz, the only signal that shows a 
Zj power above the Icr confidence level - regardless of the en¬ 
ergy band of the search - is at vaoii = 30.252122(22)Hz (see 
lower plot of Fig. O. The measured fluctuation in our tests is 
Zj -22- 39, consistent with a pulsed fraction of pf = 11.5(9)%. 
The overall signficance of the signal is low: on average, the 
probability of obtaining such a peak by chance is ~ 23%. Fur¬ 
thermore, the shift in frequency with respect to the first XMM- 
Aewfon observation is too large (v,-vaoii ~ 23.4 Hz) for the two 
signals to be associated. The implied spin-down rate would pro- 
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duce, within the usual magneto-rotational dipole-braking sce¬ 
nario in a vacuum, a magnetic held of Bdip = 6.6 x 10'^^ G. In 
particular, a held with B > 5 x lO'^ G would spin down the pul¬ 
sar at a rate >0.1 Hzyr making the detection of a very short 
period of ~ 20 ms unlikely in the hrst place. 

Regardless of the absence of an obvious candidate periodic¬ 
ity, the new XMM-Newton observation can still be used in combi¬ 
nation with the A09 data in a coherent two-dimensional Zf(v, v) 
search, i.e. allowing the test to account for the neutron star 
spin-down. In this case, the ephemeris parameters (v, v) that de¬ 
termine the phase (pj of each photon time-of-arrival, 

^ . .(tj-to? . , 

(pj = v(t - to) + V - - - ;=l,...,Aph, 

are estimated as the values that give the highest power of the 

statistic. Here, to is the event time-of-arrival counted from 
an epoch of phase zero and Nph - 10 355 is the total number 
of photons in the search, when joining the two datasets for the 
optimal choices of energy band and extraction region. 

We searched for significant Zj peaks in an interval of +5cr 
around v*, in steps of 5 x 10 ®Hz, and explored frequency 
derivatives from -1 x 10^'^Hzs^' up to a maximum of Vmax = 
-lx 10 Hz s', in steps of -2 x 10 Hz s'. The choice of 
the maximum frequency derivative is reasonable, as only three 
isolated pulsars - the Crab and the energetic pulsars PSR J0537- 
6910 and PSR 10540-6919 in the Large Magellanic Cloud, out 
of about 2400 known objects - are known to slow down at a 
rate faster than Vmax- The number of independent trials is given 
by Af = (Vmax - Vn,in)|vmax “ Vmin\T^p^J2 ~ 4 X 10^ where 
Tspan = 6.4 X 10^ s is the time span between the A09 and AOl 1 
observations. 

We found that the power of the Zj statistic increases sig¬ 
nificantly when folding the two datasets, even though the num¬ 
ber of trials is boosted by a factor of 100. The highest peak 
resulting from our search, with Z^ ~ 58, occurs at v = 
53.645805460(5)Hz and v = -2.536200(15) x 10-"Hzs-' 
{To - MID 55536.008840 days is the epoch of phase zero); 
the chance probability of the solution considering all number 
of independent trials is 0.01% (confidence level of 3.885cr). A 
total of 17 solutions, out of ~ 6 x 10^, were found at or above 
the 3cr level (with 51.5 < Z^ < 58); none was above the 4cr 
level. All high-Zj solutions have frequency derivatives within 
-9.8 X 10“" < V < -1.2 X 10“"Hzs“'. Within the scenario 
of magnetic braking in vacuum, the spin-down values imply a 
dipolar field of Bdip = (3 - 8) x lO" G for the source. 

To test the result beyond its formal significance, we repro¬ 
duced (in full resolution and over the same parameter space) 
the two-dimensional Zj (v, v) search, but folded AOl 1 events ex¬ 
tracted from a background region over the original A09 time 
series of the source. The background region was scaled to have 
roughly the same number of counts as collected for J1046, and 
two different energy bands were tested; the 0.36 - 2.2 keV of 
the original search, where the S /N ratio is optimal, and the 
3 - lOkeV energy band, which aims at excluding contamina¬ 
tion from the neutron star. The goal is to verify empirically the 
power of the Z^ statistic on a search with similar sensitivity, total 
duration, and time span of the original one, but conducted over 
the times of arrival of noise events that should not relate with the 
periodic signal at v,. 

We found that Z^ peaks of similar power are obtained on 
these tests for intervals of frequency derivative that are roughly 
of the same order as those of the original search. In the 0.36 - 


2.2 keV energy band, we found 14 solutions above the 3cr level 
(again out of ~ 6 x 10^ grid points), with Z^ ~ (51.5-61) and 
-3 X 10" < V < -8 X 10 Hz s '. The highest Zj power 

occurs at V-8 x 10“'^ Hz s“' and implies B^ip ~ 2.3 x 10" G. 

For the 3 - lOkeV energy band, 28 solutions with frequency 
derivatives in the range -7 x 10“" < v < -1.2 x 10“" Hzs^' 
were found with Zj ~ (51.5 - 60), and the highest power is for 

V-7x 10“" Hz s“', which implies Bdip ~ 7x lO" G. The tests 

show that the high Z^ power is likely to be caused by statistical 
fluctuations in a search with a large number of trials and that the 
best (v, v) solution derived for J1046 is not significant. 


4. Discussion 

4.1. Summary of results 

At present, we cannot rule out the possibility that the periodic 
signal at ~ 18.6 ms, which was only significantly detected in 
the A09 observation (pf - 12.5(1.9)%, c.l. 3cr to 4cr, depend¬ 
ing on the S IN ratio of the search; see Sect. Ira is unrelated 
to the true spin period of the neutron star. The second XMM- 
Newton observation proved inconclusive at confirming the fast 
candidate spin, providing an upper limit that is close to the lim¬ 
iting sensitivity for detecting the A09 modulation {pf < 14% at 
4cr, for P - 0.0114 - 100 s and 0.35 - 2keV). No other peri¬ 
odic signal was found to be statistically significant in extensive 
searches performed over the individual datasets (in a wide fre¬ 
quency range and for different tested energy bands). In partic¬ 
ular, pulsations with P > 0.6 s (a range relevant for AXPs, the 
M7 INSs, and RRATfl) are con strained down to pf ~ 10% (3cr, 
0.35 - 2keV: iPires etai]l2012l) . Such smooth X-ray pulsations 
are usually more frequently observed in the M7 and CCOs, as 
AXPs typically show pf > 15 - 30%. 

The coherent combination of the two EPIC-pn time series, in 
a joint periodicity search around the candidate spin and account¬ 
ing for a wide range of possible pulsar spin-down values, pro¬ 
vides a {P, P) solution with a probability of 1 in 10 000 of it be¬ 
ing spurious. Across the searched parameter space, the solutions 

above 3cr imply a dipolar magnetic held of Z?dip = (3-8)xl0"G, 

under the assumption of magnetic braking in vacuum. However, 
additional tests suggest that the timing solutions are likely to be 
below their formal statistical significance. 

In the absence of a clear period determination and good es¬ 
timate of the magnetic held, the evolutionary state and nature of 
this neutron star remain as yet to be constrained. Unfortunately, 
at the flux level of the source and with current X-ray facilities, 
pulsations are difficult to detect below the limits provided by 
present data. 

We were nonetheless able to study the spectral properties 
of J1046 with unprecedented detail and statistics (Sect. ED- 
Interestingly, we found that good fits and meaningful physi¬ 
cal parameters are only derived when the residuals at energies 
0.55 keV and 1.35 keV are taken into account by the spectral 
modelling. Independently of the exact choice of the thermal 
continuum, we found that the former residuals may result ei¬ 
ther from a local overabundance of oxygen in the Carina nebula 
(with Zo = 1.2-1.7, relative to solar) or from the inhomoge- 
nous temperature distribution on the surface of the neutron star 
(the best-fit double-temperature models have kTJ” = 30-36 eV 


' As of 2015, over 1000 sources in the ATNF pulsar catalogue (of 
which 30 are X-ray puls ars and 58 are RRATs) are known to spin at 
periods longer than 0.6 s dManchester et alJ2005h . 
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and kT^ - 120 - 190 eV). Phenomenologically, the residu¬ 
als at 0.55 keV can also be described by a Gaussian absorption 
line with EW ~ 70 - 90 eV. On the other hand, of all mod¬ 
els tested in Sect. EH the residuals at 1.35 keV are only satis¬ 
factorily accounted for by invoking a line in absorption, with 
EW ~ 120 - 400 eV (depending on the choice of the contin¬ 
uum). In this case, the fits lead to overall physically consis¬ 
tent (as well as much more tightly constrained) spectral param¬ 
eters; Nh = (2.5 - 3.3) X 10^1 cm-2, = (5 - 8) x 10^ K, 

= (0.7 - 4) X 10^^ erg s“*, and d ~ 1-3 kpc (considering the 
solutions with;^fy < 1.4). 

4.2. Comparison with other puisars 

Of all the INSs that clearly sho w thermal X-ray em ission (see, 
for example, the compilation of IVigano et akllTOl^ . J1046 dis¬ 
plays spectral properties that are quite close to those of the M7 
and CCOs - namely, the soft energy distribution, constant X-ray 
properties on a long timescale (> 15 yr), absence of significant 
magnetospheric emission, and evidence of absorption features. 
Deeper radio and y-ray limits are, however, needed to exclude a 
middle-aged rota tion-powered pulsar (RP P) or a mildly-recycled 
neutron star (e.g. lBelczvnski et al.ll2010h . 

The likely presence of J1046 in the Carina open cluster ex- 
cl udes a neutron st ar much older than ~ 10^ yr (see discussion 
in iPires et al.ll2012l) . In millisecond pulsars (MSPs), as well as 
in middle-aged and old RPPs, a considerable fraction of the X- 
ray emission arises from polar caps, which are heated by back- 
flowing charges accelerated in the pulsar magnetosphere. The 
spectral analysis shows that the X-ray emitting region in J1046 
is consistent with the surface of a canonical neutron star when 
realistic atmosphere models are adopted. This is in contrast to 
the ~ 0.1km blackbody rad ii typically derived for MSPs (e.g. 
iBogdanov et al.ll20^l201 Ih . 

Middle-aged RPPs (e.g. the Three Musketeers, with ages of 
a few ~ 10^ yr and spin-down power of a few ~ 10^“^ erg s“') are 
dominated by thermal emission and show much lower residual 
magnetospheric ac tivity (at levels withi n 0.3 % and 1.7 % of the 
source luminosity : lDe Luca et al.ll2005h . In J1046, no such weak 
hard tail (above 0.5% of the source luminosity) is present, and 
we found no compelling evidence for a second (hotter) thermal 
component, as that typically arising from heated polar caps (see 
discussion in Sect. 14.31 1. For most RPPs, the efficiency at convert- 
ing rotational energy in to X-rays is generally above rj - 10“^ 
dKargaltsev et al.ll2012h . By considering the 3cr upper limit on 
the source non-thermal X-ray luminosity, < 2 x 10^®ergs“* 
(Sect. 13.11 for the best-fit nsa model (e) and B - 10'^ G), and 
assuming a maximum age of 10® yr for the neutron star, we can 
derive a lower limit on the spin-down luminosity of the source, 
£ > 2.1 X 10^® erg s ', if the neutron star is indeed a fast rotator. 
This implies an exceptionally low efficiency ratio of < 10~® 
for J1046, if the neutron star is a typical RPP. Alternatively, if 
the signal at P, is not related with the INS spin period, the pos¬ 
sibility remains that J1046 is a fainter and more distant M7-like 
INS. 

4.3. Origin of absorption iines 


The detection of spectral features in absorption in several (non¬ 
accreting) thermally emitting neutron stars - including the 
M7, CCOs, magnetars, and Calvera - have been widely re- 


ported (e.g. ISanwal et al. 

2002; Mori et al.l 

2005t lHaberlll2007t 

iMcLaughlin et al.l 20071 

Zane et al.l 1201 1; 

Tiengo et alJl2013h. 


The physical interpretation is not unique. The features can be due 
to atomic transitions in the magnetised atmosphere or condensed 
surface of the neutron star, or be cyclotron lines generated in a 
hot ionized layer near the surface. Importantly, the central en¬ 
ergies of the features have been used to give esti mates on the 
state o f magnetisation of the neutron star. Recentlv. lVigano et al.l 
(l2014h showed that in so me cases (as for the M7 RX J0806.4- 
4123; iHaberl et al.l 2004ll . rather than an actual transition be¬ 
tween energy levels, the spectral deviation may be induced sim¬ 
ply by the inhomogeneous temperature distribution on the sur¬ 
face. 

In J1046, we verified that a double-temperature spectral 
model is only a plausible (physical) solution when we intro¬ 
duce an absorption feature at 1.35 keV. Even in this case, the 
reduced chi-square values are still high (with respect to single¬ 
temperature models) and the radiation radii are uncertain. Pos¬ 
sibly, the poor fit reflects the inadequate description of the 
anisotropic temperature distribution on the surface of a neutron 
star by the simple superposition of two atmosphere models. 

If we assume that the feature at energy e ~ 1.35 keV is the 
fundamental electron cyclo tron absorption, simil ar to the case of 
the CCO IE 1207.4-5209 (iBignami et al.ll2003h . we can derive 
an estimate of the magnetic field on the surface for J1046. Ac¬ 
cording to the relation ecyc = 1.16(B/10*'G)(l-i-Zg)^', we obtain 
Bcyc ~ 1.5 X 10" G for the gravitational redshift of a canonical 
neutron star, Zg ~ 0.3. 

4.4. Association with a runaway star in Carina 

A runaway massive star in the Carina nebula, MJ 218 (listed in 
the UCAC4 catalogue as 153-055048) was recently proposed to 
have been associated with the progenitor of J1046 in a binary 
system, which was then disru pted by the supernova explosion 
that created the neutron star (iNgoumou et al.ll2013h . At a dis¬ 
tance of 2.3 kpc, the 18' angular separation of the two objects 
on the sky implies a physical distance of ~ 12 pc. The total 
proper motion of the runaway star corresponds to a t ransverse 
speed of Vt = 95 + 52kms“’ (IZacharias et al.ll2004t) . Assum- 
ing typical masses of 1.4 Mo and 9.6Mo (spectral type B1.5 V; 
iPecaut & Mamaie3l2013h . respectively, for the neutron and run¬ 
away stars, the estimated flight time is of fyn ~ (1.1-3)x lO'^yr, 
which is roughly consistent with the absence of a supernova 
remnant. Although uncertain, the implied velocity of J1046, 
Vt ~ 650 + 360kms“', if previously associated with MJ 218, 
is also consistent with the distributi on of transverse sp eeds ob¬ 
served in non-recycled radio pulsars (iHobbs et al.ll2005h . 

4.5. Thermai evoiution 

Based on our best-fit spectral solutions, we can investigate the 
thermal state of the neutron star amidst the several INS popu¬ 
lations for the first time. Eigure |5] shows a representative sam¬ 
ple of INSs, in a cooling-age diagram. The temperature £“ is 
the apparent effective surface temperature, as derived from ei¬ 
ther atmosphere or blackbody models. Whenever possible, age 
estimates and intervals are derived from proper motion mea¬ 
surements (flight time) or associations with supernova remnants; 
otherwise the pulsar spin-down age, uncertain by a factor of 
two, is adopted. The data on RPPs (14 objects including the 
Crab and Vela pulsars, as well as the Three Mus keteers and the 
old pu lsar PSR J2043H-2740) are collected from I Yakovlev et El 
(l2008b . The five CCOs in the diagram are the neutron stars in 
the Cas A, Puppis A, G296.5, Kes 79, and HESS J1731-347 su- 
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Fig. 5. Cooling-age diagram for different groups of isolated neutron 
stars (data points: see legend). Theoretical cooling curves for a canon¬ 
ical neutron star are shown as lines (see text). No effects of magnetic 
field decay are considered. 


pernova remnants (iBogdanovI 120141 : iKlochkov et al.llT015L and 
references therein); in particular, their emission can be inter¬ 
preted as radiation from the entire surface, using either hydrogen 
or carbon atmosphere models. The locus of magnetars (SGRs 
and AXPs) is represen ted by the seven objects investigated by 
iKaminker et al.l (120091 and references ther ein). Data points for 
the M7 are according to the compilations of iKaplan et al.l (1201 Ih 
and IVigano et al.l (l2013h . For these sources, characteristic ages 
are shown as upper limits on true ages in the absence of a kine¬ 
matic estimate. 

Following our previous discussion, we adopt J1046’s flight 
time as a lower limit for its age estimate; a conservative value 
of 10^ yr is taken as an upper limit, based on an association with 
the Carina nebula. The temperature range is as reported in Sec¬ 
tion 14.31 For Calvera, w e adopt the tempe rature of the best-ht 
atmosphere (nsa) model (IZane et al.ll201 Ih and th e INS charac¬ 
teristi c age (again uncertain by a factor of two, iHalpern et al.l 
120131) . 


The set of theoretical cooling curves in Figure |5] corre¬ 
spond to four representative f a milies , which are descr i bed in , 
for example, ^^isskopf et al.l (1201 lb ; lYakovlev et al.l (1201 Ih : 
iKlochkov et al.l (12015111 In a nutshell, the so-called standard 
neutrino candle, represented in Figure |5] by a solid line, shows 
the thermal evolution of a canonical neutron star (of mass M - 
1.5 Mq and radius R - 12 km, with a nucleon core), where cool¬ 
ing is governed by the modified Urea process, and an iron en¬ 
velope is assumed (no enhanced cooling in the hrst 100 years). 
The dotted line shows the same star affected by strong proton 
superfluidity in the core, which highly suppresses the Urea reac- 


^ As noted bv iKlochkov et alj ll201.5h . these curves are mostly aimed at 
explaining the residual heating observed in CCOs. Colder neutron stars, 
like some RPPs in Figure [5] may be regarded as sufficiently massive 
stars with higher neutrino luminosity owing to the onset of the direct 
Urea process in the core. 


tions and results in a hotter neutron star at f > 100 yr. The dashed 
and dot-dashed lines are the corresponding (non-superfluid and 
superfluid) models for a star with a full amount of accreted 
light elements (in particular, AM ~ IO^^Mq of carbon) in the 
blanketing envelope. The presence of light elements in the en- 
velope of the neutro n star increases the thermal conductivity 
dPotekhin et al.l 119971) : as a result, the light-element envelope 
regulates and slows down the cooling in the first ~ 10"^ - 10^ yr, 
depending on the amount and composition of the accreted mate¬ 
rial. The effect of the magnetic field on cooling is not taken into 
account, which is a valid approximation for B < lO'^ G. 

The strong thermal emission (high temperature at a given 
age) observed in magnetars and the M7, with respect to the 
standard cooling theory and other INS groups, cannot be ex¬ 
plained without taking into account the effects of magnetic 
field decay as an additional source of heating o f the neutron 
star c rust. Magneto-thermal evolutionary models dVigano et al.l 
1201 3l) show that, for strongly magnetised INSs, field dissipa¬ 
tion keeps the stellar crust hot for a longer time than expected 
from standard cooling (see. I Aguilera et aDl2008l and references 
therein); moreover, strong fields at birth (B > 10'"^G) can sig¬ 
nificantly brake the neutron star spin to an asymptotic value in 
a relatively short timescale (~ 10^ yr). Interestingly, whereas for 
the bulk of the neutron star population such effects are negligi¬ 
ble, the model implies an evolutionary link between the M7 and 
magnetars. 

On the other hand, consistent with the standard theory, CCOs 
are young neutron stars that show the slow cooling in the 
neutrino-cooling era typical of sources with light-element ac¬ 
creted envelopes. For these neutron stars, the turning point to¬ 
wards the steep drop in temperature of the photon-cooling stage 
takes place earlier than for neutron stars with an iron envelope. 
As a result, old CCOs may be invisible to X-ray observations 
after ~ 1 Myr, where luminosities may have dropped below 
~ 10^' ergs ', and surface temperatures are colder than ~ 20eV 
to 30 eV. Nuclear burning in the envelope, and hence a change of 
composition, might nonetheless slow down t he cooling in the l ast 
stages of the thermal life of the neutron star dPage et alj|201 ll) . 

Overall, the cooling status of Calvera appears consistent with 
a scenario where this neutron star evolved from a CCO (see 
Sect. [B. However, Calvera’s real surface temperature may be 
considerably softer than the adopted value in Figure |5] kT ~ 
98 eV, as the source shows a clear indication o f a double¬ 
temperature energy distribution (IZane et alJl201 ll) . In particu¬ 
lar, the fit with a double thermal (blackbody or nsa) model con¬ 
strains the apparent temperature of the cooler component within 
r“ ~ (6 - 9) X 10^ K, which is very similar to the interval es¬ 
timated for J1046. Given its resemblance to Calvera, and if the 
fast spin and indication of a low magnetic field in J1046 are to be 
confirmed, the INS in Carina may be another potential candidate 
for the elusive class of evolved anti-magnetars. 


5. Summary and conclusions 

Following the intriguing detection of a short periodic signal i n 
the neutron star 2XMM J104608.7-594306 (iPires et al.ir2^ , 
we targeted the source again with XMM-Newton, aiming at better 
characterising its spectral energy distribution and timing prop¬ 
erties. The new XMM-Newton dataset unfortunately proved in¬ 
conclusive in confirming the fast spin of the neutron star. The 
derived 4cr upper limit on the pulsed fraction, pf < 14%, shows 
that the observation is just at the limiting sensitivity to detect the 
modulation found previously. 
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Without a clear determination of the spin period and a robust 
estimate of the magnetic field of the neutron star, its exact na¬ 
ture remains open to interpretation. The INS 2XMM J104608.7- 
594306 may be similar to Calvera, a neutron star for which 
the scenario of an evolved anti-magnetar, within the frame¬ 
work of the magnetic field burial by fallback of supernova 
debris, has been discussed. The overall spectral properties of 
2XMM J104608.7-594306 and its likely presence in the Carina 
open cluster disfavour a recycled object or a standard evolution¬ 
ary path, typical of that of a rotation-powered pulsar. However, 
better age estimates (e.g. through kinematic studies), as well as 
deeper radio and y-ray limits, are required to further constrain 
the evolutionary status of this neutron star. In particular, stud¬ 
ies in the optical and near-infrared may offer the opportunity to 
assess properties which, for J1046, are not assessible by other 
means. 

The prospect of hnding more thermally emitting isolated 
neutron st ars with the upcom i ng all-sky survey X -ray mission, 
eROSITA (iPredehl et al.ll2014l iMerloni et alJl201^ . is of course 
of much interest. eROSITA will scan the X-ray sky at unprece¬ 
dented flux levels, making it possible to recognise orphaned 
CCOs and M7-like neutron stars through their residual thermal 
emission, either among the known radio pulsar population or 
even if these elusive objects are intrinsically radio silent. 
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